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ABSTRACT 

In this paper, we present an original observational approach, which combines, for the first time, traditional 
speckle imaging with image post-processing to obtain in the optical domain diffraction-limited images with high 
contrast (10 -5 ) within 0.5 to 2 arcseconds around a bright star. The post-processing step is based on wavelet 
filtering an has analogy with edge enhancement and high-pass filtering. Our I-band on-sky results with the 
2.5-m Nordic Telescope (NOT) and the lucky imaging instrument FASTCAM show that we are able to detect 
L-type brown dwarf companions around a solar- type star with a contrast AI~12 at 2" and with no use of any 
coronographic capability which greatly simplifies the instrumental and hardware approach. This object has 
been detected from the ground in J and H bands so far only with AO-assisted 8-10 m class telescopes (Gemini, 
Keck), although more recently detected with small-class telescopes in the K band. Discussing the advantage and 
disadvantage of the optical regime for the detection of faint intrinsic fluxes close to bright stars, we develop some 
perspectives for other fields, including the study of dense cores in globular clusters. To the best of our knowledge 
this is the first time that high contrast considerations are included in optical speckle imaging approach. 
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1. INTRODUCTION 

To date, a high number of ambitious scientific cases of modern galactic and extragalactic astronomy rely on the 
need of very deep and sharp images able to unveil faint details. Such goals stress the importance of observations 
with high sensitivity high dynamic range and high spatial resolution. In the last ten to twenty years, these 
requirements have been fulfilled with the development of a new class of large 8-10 m telescopes equipped with 
adaptive-optics systems to compensate for the atmospheric turbulence that strongly degrades the observational 
needs mentioned above. Space-based observatories such as the 2.5-m HST - or the future JWST - have been so 
far the only alternative to AO-based instruments to obtain highly detailed images. 

In the field of low- mass stars, brown dwarfs and exoplanets, the potential of this approach is well recognized and 
extremely vast: direct imaging at high spatial resolution and high contrast permits to reveal multiple systems 
and planet candidates with longer periods (~0.1") than what typically probed with radial velocities, to char- 
acterize the circumstellar environments of young forming planetary systems, 1 or to determine unambiguously 
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dynamical masses by following over few years the orbit of close brown dwarfs binaries in order to calibrate the 
mass-luminosity relationship for substellar objects. 2 ' 3 In other words, the determination of masses of very 
low-mass and planetary mass objects, as well as the establishment of multiplicity statistics and orbital properties 
are among the most fundamental parameters that can be measured experimentally using the "direct imaging" 
approach, which in return will help us to constrain star and planet formation theories. 

From the ground, high angular resolution (HAR) observations rely on adaptive-optics systems, which are pri- 
marily designed to operate in the near-infrared regime (J, H, K bands). In the optical regime (A<1 /im), com- 
plementary data to the infrared ones are necessary for a full characterization of the spectral energy distribution, 
key to the determination of effective temperatures and bolometric luminosities. However, one is generally limited 
at these wavelengths by the poor correction of conventional near-IR optimized AO systems, unable to deliver 
diffraction limited resolution over significant (>5 // ) field of view. 

Optical speckle imaging techniques 4, 5 share similar goals with near-IR AO-based instruments. The strong advent 
of fast readout and low readout-noise CCD has favoured the development of "Lucky Imaging" instruments able 
to deliver optical diffraction-limited images in which faint companions and close binaries can be detected at vis- 
ible wavelengths, offering complementary scientific outcomes to near-IR observations. In this paper, we present 
new results obtained with the FASTCAM instrument 6 developed at the Instituto de Astrofisica de Canarias 
which address both high angular resolution and high dynamic range issues. After a short presentation of the 
instrument, we give an overview of its scientific capabilities, stressing on the potential for high dynamic range 
imaging in association with image post-processing. 

2. INSTRUMENT DESCRIPTION 



2.1. Exploiting diffraction-limited speckles 

As soon as the size of the telescope overcomes the Fried parameter ro that describes the atmospheric turbulence 
at a given wavelength, a point-like source image will appear as a rapidly spatially and temporally varying speckle 
pattern within the characteristic timescale of the turbulence. In long-exposure images, this effect will smear 
out the diffraction-limited point spread function (PSF), leading to seeing-limited images for which the spatial 
resolution becomes independent of the size of the telescope (see Fig. 1). 

Speckle imaging has been one of the pioneering high resolution imaging technique able to take advantage from 
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Figure 1. Left: short-exposure image of the binary star GN Tau. The pattern is composed of different speckles with 
different brightness. Center: The seeing limited image of the source after a long exposure. The various speckles add 
randomly, resulting in an diameter independent FWHM of the PSF. Right: The shift-and-add plus frame selection image. 



the high spatial resolution information contained in speckles, each of them representing an image of the source 
at the diffraction limit of the telescope. Speckle imaging is not a new technique: it has been implemented since 
the seventies by different groups, 7,8 but it is only in the last years that optical speckle imaging made real 
advance in astronomy (see McCay et al., these proceedings) essentially thanks to the very strong improvement 
of optical CCD detectors. Indeed, because this technique relies on the acquisition of short-exposure images, 
the readout-noise of the detector becomes the dominant noise source limitation. Going towards fainter sources 
requires to decrease as much as possible the readout noise, which becomes now possible with the new generation 
of fast readout CCD detectors. 

In brief, speckle imaging stores separated short-exposure images lasting few tens to few hundred of milliseconds, 
each of them presenting a random realization of the atmospheric turbulence. The next step is to recenter all 
the images of the dataset along the brightest speckle - determined by some metrics, e.g. the brightest pixel in 
a given region - using a shift- and- add technique to collapse the image cubes. The final PSF is composed of a 
diffraction-limited core resulting from the co-addition of the brightest speckles in each frame, surrounded by a 
halo - the "seeing disk" - caused by the superposition of the other speckles in the image, and similarly present 
in classical AO images due to imperfect correction of the turbulent wavefronts (see Fig. 1). 
By applying a selection of a given percent of best frames in the sequence, it is possible to improve at the same 
time the spatial resolution (decreasing the FWHM of the PSF) and reduce the relative contribution of the halo. 
This is because the best frames in the sequence corresponds to the lowest turbulence effect of the atmosphere 
(i.e. larger ro). This is illustrated in the right side of Fig. 1. "Best frames" is intended here as those frames 
where most of the intensity in the image is concentrated in one given speckle. 

2.2. The optical speckle camera FASTCAM 

FASTCAM is a collaborative project developed by different teams at the Instituto de Astrofisica de Canarias, 
Canary Islands and at the Universidad Politecnica de Cartagena, in Spain. The current view of the instrument 
after the opto-mechanical design phase is shown in Fig. 3. The FASTCAM instrument has a great versatility and 
can equip different telescopes of the ORM observatory* such as the 2.5-m Nordic Telescope (NOT) t or the 4.2-m 
William Herschel Telescope (WHT) Most of the time, FASTCAM is installed as a common- user instrument at 
the 1.5-m Carlos Sanchez telescope (TCS). 



*ORM: Observatory of Roque de los Muchachos 
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Figure 2. View of the FASTCAM camera on its mechanical support in the standard common-user instrument con- 
figuration. The speckle camera is seen on the right part of the image. The mechanical support holds also two filter 
wheels. 



Our speckle camera is based on a commercial fast readout (EM) CCD detector developed by Andor Technol- 
ogy. Basically, this new generation of optical detectors present the strong advantage of performing the "Electron 
multiplication" (EM) of the signal before the readout of the detector, which means that, proportionally, the 
readout noise can be significantly reduced in each pixel compared to useful signal, down to ~le- rms. The cur- 
rent detector is 512x512 pixels, with possibility of windowing the array to speed up the readout. The individual 
integration time can be effectively set down to 30 ms. The quantum efficiency is quite high, peaking at 95% at 
0.55 fim. 

The optical system of FASTCAM is based on relatively simple doublet lenses that permit to deliver various 
pixel scales, namely 70 to 50mas/pix at the TCS, 30mas/pix at the NOT and 13.5mas/pix at the WHT. Cur- 
rently, no atmospheric dispersion compensator (ADC) is used with FASTCAM, although a dedicated ADC is 
being designed. However, this issue is not very constraining as long as we observe as low airmasses or in narrow 
band filters. A motorized filter wheel permits to install the majority of the wide band and narrow band optical 
filters typically used for imaging at visible wavelengths. 




Figure 3. Screen shot of the computer interface controlling Fastcam (left) and the data reduction software (right) 
developed by the University of Cartagena, Spain. 



FASTCAM benefits from its own software interface to control the instrument at the telescope, as well as from 
a dedicated data reduction software used for both an accurate data reduction of the speckle images and for a 
quick view of the data at the telescope. The control software and data reduction pipeline were developed at the 
University of Cartagena. 



3. SCIENCE WITH FASTCAM 
3.1. Multiple systems at high angular resolution 

Speckle imaging has been primarily used for imaging close sub-arcsecond multiple systems from the ground in the 
optical, while they were generally accessible only using AO systems in the near-infrared. Monitoring of binary 
stars is one of the main objective of FASTCAM, in particular those with small angular separation and relatively 
close in order to determine in a reasonable time scale the orbital parameters of the binary, which ultimately 
provide dynamical masses of the system. 

Different observing campaigns have been conducted at the three different telescopes, reaching close to diffraction- 
limited images in each occasion. Fig. 4 illustrates some of the results obtained with FASTCAM in R and / bands 
at the 4.2-m WHT. For instance, the binary system CHR 36 with a separation of 0.063" was resolved in the /-band 
filter centered at 790 nm using integration times of the order of 50 ms. 

Speckle imaging is not only limited to a small field-of-view in which to observe binaries. High resolution imaging 
can be performed on a larger field-of-view, which can be illustrated with the notion of isoplanatic angle. The 
frame selection process keeps those "best" frames corresponding to the lowest atmospheric turbulence (i.e. larger 
ro) for which the isoplanatic angle is larger than the average one obtained when 100% of the frames are main- 
tained. Hence, the finest frame selection can guarantee isoplanatic angles of 30" or larger, significantly larger 
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Figure 4. Close binaries images in the I-band and R-band using FASTCAM at the 4.2-m William Herschel Telescope. 

than the average isoplanatic angle of in normal seeing conditions. This can be successfully exploited for 
imaging crowded fields in dense cores. 

Fig. 5 compares two images of the triple system UZ Tau with comparable resolution. In the center part, UZ 
Tau is imaged at 2.16 /im with the AO-equipped NACO instrument at the VLT. 9 On the right is shown the 
optical image counterpart obtained at 0.8 /am with FASTCAM at the 2.5-m NOT. The relative brightness of 
each component at 0.8 //m is similar to the infrared case, which provides indications on the spectral type of the 
system components. The AB separation is ~3.5" while the BC separation is ^0.36". 
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Figure 5. Comparison of infrared and optical images of the triple system UZ Tau obtained with NACO (left) and 
FASTCAM at the NOT (right). 



3.2. The question of high dynamic range at short angular separations 

On the right side of Fig. 5 is shown a close view of the resolved system UZTauBC, which at the same time 
highlights a critical issue of lucky imaging: around each of the two components the seeing halo is likely masking 
any faint structure or companion that may lie at short angular separations. For classical high contrast long- 
exposure imaging, this question is traditionally solved with the use of a coronograph placed in the intermediate 
image plane that strongly reduces the flux contribution of the star on-axis and consequently of the halo residual 
left by the AO. In absence of any coronograph, a possible approach with speckle imaging is to reach as much 
as possible the true diffraction-limited case that corresponds to the highest Strehl ratios. This can be done 
with speckle image reconstruction for which several methods have been developed in the past, such as bi-speckle 
inter ferometry, speckle holography, blind deconvolution etc... 

A second approach, that we decided to adopt in this study, is based on numerical post-processing of the lucky 
image in order to filter out low spatial frequencies and enhance other ones, in particular those related to faint 
nearby companions. This approach is based on the wavelet decomposition of the image - a standard and 



routinely used image processing technique in computer sciences - that permits to identify the best trade-off 
between the optimal spatial frequency at which the companion is revealed and the accessible signal-to- noise. 
Here, we implemented an iterative wavelet filtering described by the relationship 1=1-1 *k, where k is the size 
variable wavelet kernel. The wavelet filtering process is illustrated in the sequence of Fig. 6. The first image shows 
the FASTCAM/NOT image at 0.8 jam of the triple system Haro6-37. The A component is actually a double 
system with a fainter companion Ab at 0.36" from Aa. In the direct lucky image, this companion can be hardly 
detected in the PSF halo surrounding Aa. The four successive images show the multi-resolution decomposition 
of the first image, with the kernel size changing by a factor 2 from the second to the fifth image. By comparing 
the spatial scale at which Ab is detected and the visual SNR of the detection, it appears that the second image 
of the sequence is the optimal trade-off between those two factors, while in the two last cases we can see that 
the spatial scale at which the image is filtered results in the partial or complete removal of the companion. The 
rebound dark zone surrounding the various components is an intrinsic feature of the filtering process. 
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Figure 6. Wavelet decomposition of the lucky image of the system Haro 6-37 observed with the NOT. See text for details. 



In practice, the wavelet decomposition is substituted with wavelet filtering in which only the optimal component 
of the decomposition is maintained as the final post-processed image. In such a way, we benefit from a relatively 
simple filtering algorithm to enhance the detection of faint and nearby point-like sources. The effect of the wavelet 
post-processing on faint extended source is not treated here, although it certainly present a strong interest. 

3.3. The case of HD 130948 system 

In order to evaluate the potential of this approach, we applied the method presented above to the case of 
HD 130948 in order to obtain a high contrast image of the system. HD 130948 is a main sequence G star at ~20 pc 
hosting a brown dwarf binary at a separation of ~2.5 // from the primary. HD 130948 BC is a benchmark for the 
independent determination of dynamical masses of sub-stellar objects, which is a key step for the calibration of 
the corresponding mass-luminosity relationship. The binary was initially discovered by Potter et al. (2002) 10 in 
the near-infrared using the 8-m GEMINI telescope coupled to adaptive optics, and consequently reported as a 
pair of L2±2 sub-stellar objects. 
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Figure 7. Left: Original detection with the 8-m GEMINI telescope; Center: Direct lucky image with the 2.5-m NOT; 
Right: HD 130948 system after wavelet image post-processing. In all the images, North is up and East is left. - Adapted 
from Labadie et al. (2010) 11 . 



HD 130948 was observed with FASTCAM on the 2.5-m Nordic Telescope at two different epochs, in May 2008 
and July 2008. We acquired a total number of, respectively, 5xl0 4 and 10 5 images in the / filter. The night 
conditions were good and excellent for the May and July runs. Fig. 7 shows respectively the original detection, 
the direct lucky image and the result after post-processing for the July run. In the last image, we clearly 
detect a point-like source east from the primary with SNR^IO, which position is consistent with that none for 
HD 130948 BC. In the NOT image, HD 130948 BC is not detected as a binary, but as a point source. From a 
posterior detailed study of the brown dwarf binary orbit 12 , the binary separation at the time of observation was 
^35 mas, beyond the actual diffraction limit of the telescope necessary to resolve the pair. 

Different tests have been conducted in order to confirm the detection and avoid a false-detection of a speckle. 
A first test has consisted in splitting the July data in three different time series in order to check if the brown 
dwarf pair remains detectable at the same location, oppositely to time varying speckle. This is illustrated in the 
sequence of Fig. 8, which presents observations at three different moments in an annulus of appropriate size. At 
the position of the white arrow is shown the location of HD 130948 BC, which is detected in all the three time 
series, while other bright features are not present in all the three sequences. This confirms that the detected 
point-source is real. In the last image, the detection of HD 130948 BC has less signal-to-noise than in the other 
two cases: this is because the frames selection for the last time serie resulted in a poorer Strehl ratio - probably 
because of degrading weather conditions - which has spread out the flux of the brown dwarf. 






Figure 8. Three different time series of the HD 130948 BC observations from July 2008, limited to , 
white arrow indicate the location of the brown dwarf in the image. 



2.5" annulus. The 



The ultimate test for the detection of the companion consists in the second epoch observation of May 2008 taken 
with FASTCAM. The same data reduction procedure is adopted. Despite the lower number of available frames, 
the detection of HD 130948 BC is also confirmed in the May dataset, although the PSF is slightly more spread 
out than in the July data (cf. Fig. 9). The two images have been normalized to the peak intensity and can be 
directly compared. 
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Figure 9. Second epoch observation of HD 130948 BC 



In order to evaluate the impact of the wavelet filtering on the photometry and astrometry of the object, we have 
simulated several companions at a known distance from the primary and with a known magnitude difference 
(e.g. Fig. 10). After having convolved this model with the direct lucky image and applied the filtering process, 
we measure the position of the photocenter and the peak intensity of the component with respect to the input 
model. We could verify that the astrometric position was not altered by more than ~0.2 pixels by the filtering 
process, and that the photometric variation is negligible down to contrasts of ~10 -5 . This result was expected 
theoretically since the wavelet filtering process is a linear one, which normally keeps the contrast value in the 
different regions of the image. 




Figure 10. Astrometric and photometric check of the effect of the wavelet filtering process 

This work on HD 130948 BC has permitted us to derive the first measurement of the combined photometry in 
the /-band and examine the astrometric parameters that reveal orbital motion. 11 These results suggest that 
optical speckle imaging has sufficient potential for high dynamic range imaging, even with a 2.5-m telescope, to 
study faint and close sub-stellar companions in a wavelength regime not efficiently accessible with the current 
AO systems. 

3.4. Exploring dense cores in globular clusters 

Another strong application of speckle imaging coupled to image post-processing is high resolution and high con- 
trast imaging in dense cores of globular clusters. These regions are composed of compact groups of thousands of 
stars in a wide range of magnitudes. The deepest and sharpest observations of these objects are achieved with 
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Figure 11. HST observations of the M15 cluster and its core at 700 and 800 nm. See text for details. 



space-based observatories like HST able to couple high sensitivity and high spatial resolution in the optical. 
In Fig. 11, we show two HST optical images of the M 15 globular cluster. The image on the left side is a wide 
field-of-view image of the cluster taken in the I band (100 mas pixel scale) where the dense core is identified in a 
15" x 15" region. The image on the right is a zoom on this region taken at shorter wavelength and 50 mas/pix with 
HST/STIS. In the center of such a region, the aggregation of an elevate number of stars results in a diffuse halo 
of light able to mask the fainter stars, especially in the red part of the optical spectrum. Although FASTCAM 
is not comparable to HST in terms of sensitivity, the high spatial resolution and high contrast capabilities of 
out instrument help us to detect from the ground a very large fraction of objects in the I band that have optical 
counterparts detected in B and V bands with HST 13 . 

In Fig. 12 are shown two FASTCAM views of the M15 core, with the same field-of-view as in the right im- 
age in Fig. 11. The left view shows the direct lucky image with 7% frame selection (for an effective integration 
time of ~7s; 0.03^/pixel) taken with the 2.5-m Nordic Telescope. A coarse and direct comparison with the HST 
image shows the sharpness of the FASTCAM point-spread function, which permits to detect relatively close and 
faint stars. In the right view of Fig. 12, the strong halo-suppression effect of the wavelet filtering applied to the 
M15 image reveals a high number of faint point sources with contrasts down to ~5x 10 -4 - 10 -3 compared to the 
brightest stars in the field. The completeness magnitude in the I-band is ~19. The number of sources detected 
at 3-<j with FASTCAM is above 1000 in a radius of 8" from the centre of M15. These source were previously 
known from HST images at shorter wavelengths. 




Figure 12. FASTCAM observations of the M15 core and corresponding wavelet filtering result. See text for details. 



4. CONCLUSIONS AND IMMEDIATE PERSPECTIVES 

We have presented in this paper first observational results with the speckle imager FASTCAM in the context of 
high angular resolution and high dynamic range imaging with the speckle imaging technique from the ground, in 
a spectral domain where classic AO is not efficient. We are able to routinely obtain close to diffraction-limited 
optical images {I band) both with the 2.5-m Nordic Telescope and the 4.2-m William Herschel Telescope on the 
Observatory of Roque de los Muchachos, in the Canary Islands. 

Our results demonstrate high sensitivity in crowded fields (7^19), while high contrast of ~12 mag is obtained 



in the vicinity (~l-2") of bright sources, with a ~2-3 mag gain due to post-processing wavelet filtering of the 
images. This permitted to detect, without resolving so far, the brown dwarf binary HD 130948 BC considered as 
a benchmark for the independent determination of dynamical masses. The property of unsaturated PSF in the 
FASTCAM images permits us, as opposed to classical long-exposure imaging, to carry high accuracy astrometric 
measurement of brown dwarf binaries and high accuracy relative photometry (see Femenia et al., SPIE Adaptive 
Optics Systems II (Conference 7736) proceedings). 

As the wavelet post-processing filtering can presumably affect the detection of extended sources, we plan in 
a near future to explore the option of more classical image reconstruction with speckle imaging, based on speckle 
holography or phase diversity techniques, in order to improve the dynamic range in the image without filtering 
out precious spatial frequencies. Ideally, these techniques could be applied to the case of imaging of edge-on 
protoplanetary disks in scattered light at optical wavelenghts, for which nearby (~5") stars could be used as 
guide stars during the observations and PSF references for image reconstruction purposes. 
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